Short-time dynamics of random-bond Potts ferromagnet 
with continuous self-dual quenched disorders* 



Z.Q. Pan,! jjp Yingi',2,| and D.W. Gu^ 

'Physics Department, Zhejiang University, Hangzhou 310027, P.R. China 
^ Universitat Siegen, D-57068 Siegen, Germany 
^Hangzhou Science and Technology University for Adult, Hangzhou 310012, PR China 



o 
o 

(N 



o 



o 



> 
o 
cn 

m 
o 



o 



We present Monte Carlo simulation results of random-bond Potts ferromagnet with the Olson- 
Young self-dual distribution of quenched disorders in two-dimensions. By exploring the short-time 
scaling dynamics, we find universal power-law critical behavior of the magnetization and Binder 
cumulant at the critical point, and thus obtain estimates of the dynamic exponent z and magnetic 
exponent rj, as well as the exponent 9. Our special attention is paid to the dynamic process for the 
g=8 Potts model. 

PACS numbers; 75.40.Mg, 64.60.Fr, 64.60.Ht 

Keywords: Short-time dynamics, Potts model. Critical exponents. Universality 
Keywords: 



I. INTRODUCTION 

Ferromagnetic systems with quenched randomness are 
studied intensively for their critical properties affected by 
an addition of disorders ||, |j . The Harris criterion was 
claimed 27 years ago that if the specific heat critical expo- 
nent a of the pure system is positive, then quenched dis- 
order is a relevant perturbation at the second-order crit- 
ical point and it causes changes in critical exponents 
Then, following the work of Imry and Wortis ^ who ar- 
gued that quenched disorder could smooth of a first- order 
phase transition and thus produce a new second-order 
phase transition, the introduction of randomness to pure 
systems originally undergoing a first-order transition has 
been comprehensively considered @, H). The theory 
was firstly numerically checked with the Monte Carlo 
(MC) method by Chen, Ferrenberg and Landau (CFL) 
PI who studied the 8-state Potts model with a self-dual 
random-bond disorder (see Eqs.(||) and (^). Recently 
most of the work for such disordered systems is carried 
out by intensive MC simulations [|o[ y, ^, 
for examining how a phase transition is modified by the 
quenched disorder coupling to the local energy density 
for the g-state random-bond Potts ferromagnet (RBPF) . 

The g-state RBPF, an interesting framework to study 
the influence of impurity on pure systems, is described 
by the Hamiltonian, 



f3H 



E 

<i,3> 



K,, > 0, 



(1) 



where f3 — l/ksT is the inverse temperature, the spin 
a can take the values \,---,q, 5 stands for the Kro- 
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neckcr delta function, and the nearest-neighbor interac- 
tions < i,i > are considered. Usually the dimensionless 
couplings Kij are selected from two positive (ferromag- 
netic) values Ki and K2 = rKi, with a strong to weak 
coupling ratio r = K2/K1, called a disorder amplitude, 
according to a bimodal distribution, 



P{K) = p5{K - Ki) + (1 - p)5{K - K2) 



(2) 



With p — 0.5, the system is self-dual and the exact criti- 
cal point can be determined by Pq|, 



(e^^ - l)(e 



K' 



1) 



(3) 



where Kc and K'^^ are the corresponding critical values of 
Ki and K2, respectively. The value r = 1 corresponds 
to the pure case where the critical point is located at 
Kc = log(l -I- and the phase transitions are first- 
order for q > A. As a > for all g > 2 the randomness 
acts as a relevant perturbation and, for g > 4, it even 
changes the nature of the transition from first to second 
order. 

Up to present, the MC studies are usually carried out 
in equilibrium and concentrated on disorder amplitudes 
in the range of r = 2-20 adapted to numerical analy- 
sis. They give good estimates of the disordered fixed 
point exponents |2^ . In Table 1 we present the mag- 
netic scaling index 77 = 2/3/i/ of the 2D 8-state RBPF ob- 
tained by different groups ^ |l^, |o[ |2^ . In a recent 
work, Olson and Young (OY) |l5| used their special self- 
dual disorder distribution and performed the MC study 
of multiscaling properties of the correlation functions for 
several values of q. Their results for the magnetic ex- 
ponent are very interesting to examine the universality 
of the RBPF in the crossover regime from the pure fixed 
point to a percolation-like limit [po| . Cardy and Jacobsen 
p7t studied the RBPF based on the connectivity transfer 
matrix (TM) formalism, and their estimates of the criti- 
cal exponents lead to a continuous variation of 77 with q, 
which is in sharp disagreement with CFL's results M. 
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TABLE I: Magnetic exponent r; estimated by different groups 
for the 2D 8-state RBPF. 



Authors 






r 


V 


Technique 


CFL|9| 






2,10 


0.236(4) 


MC 


Olson and Young 


15 


Px{x) 


0.319(7) 


MC 


Cardy and Jacobsen 


i 


2 


0.284(4) 


TM 


Chatelain and Berche 


18 


1 10 


0.304(6) 


MC 


Chatelain and Berche 


19 


1 10 


0.301(1) 


TM 


Picco[2C] 






10 


0.306(2) 


MC 


Ying and Harada 


21 


1 


10 


0.302(5) 


STD 



In this paper, we present a MC study to verify the 
dynamic scahng of the RBPF and estimate the critical 
exponents in the short-time dynamics (STD) |2^, for 
the OY's continuous distribution of the disorder inter- 
actions. It is a self-dual scheme to introduce quenched 
randomness and defined by the following probability 



Pxix) 



2V9 



7r(l 



(4) 



for the Boltzmann factor x = . It is obvious that the 
model has strong disorder as the distribution function is 
nonzero everywhere at a; £ (0, 1) {K —^ oo, K —^ 0) and 
has a finite weight at its two limit points a; = and x ~ I. 
Another advantage of this distribution is that it is easy 
to generate random numbers with the probability Px (x) 
through 



1 



1 + ^tan(7ri?^/2)' 



(5) 



where Rx is a random number with uniform distribution 
between and 1. 

In particular, we will investigate the critical behav- 
ior affected by introduction of such continuous quenched 
randomness to clarify the universality class of the RBPF 
in the crossover regime. In simulations we mostly choose 
the model with q = 8 which is known to have a strong 
first-order phase transition without disorders. We hope 
to confirm the new second-order phase transition induced 
by the random disorder in order to show the influence of 
quenched impurities on the first-order systems and to 
check whether an Ising-like universality would be satis- 
fied from the STD approach. 



II. SHORT-TIME DYNAMICS 

For a long time, it was believed that universality and 
scaling relations can be found only in equilibrium or in 
the long-time regime. In Refs. [ p2| |2^ , however it was 
discovered that for a 0{N) vector magnetic system in 
states with a very high temperature T ^ Tc, when it is 



suddenly quenched to the critical temperature Tc, and 
evolves according to a dynamics of model A, a univer- 
sal dynamic scaling behavior emerges already within the 
short-time regime, 

Af^*^) (i, r, L, mo) = b''"^/'' M^''^ {b-'t, b^'^r, b'^L, 6^°mo), 

(6) 

where M^''^ is the fcth moment of magnetization, t,T ^ 
[T — T^jTc^ L and b are time, reduced temperature, lat- 
tice size and scaling factor, respectively. /3 and v are 
the well known static critical exponents. The xq, a new 
independent exponent, is the scaling dimension of initial 
magnetization mg. This dynamic scaling form is gener- 
alized from finite size scaling in the equilibrium stages. 
The dynamic MC simulations have been successfully per- 
formed in such non-equilibrium processes to estimate 
the critical ^oint and the critical exponents (K ,fi ^Q) 

The MC simulation begins with a study of the time 
evolution in the initial stage of the dynamic relaxation 
starting at very high temperature and small magnetiza- 
tion (mo ^ 0). For a sufhciently large lattice (L cxd), 
from Eq.(||) by setting r = and b = t^l^ , it is easy to 
derive that 



M('=)(i,mo) = i-'=''/''^M('=)(l,i^"/^mo). 



(7) 



When /c = 1 we get the most important scaling relation 
on which our measurements of the critical exponent Q are 
based, 



M{t) - mot^ 



(xo - l3/v)/z. 



(8) 



In most cases the magnetization undergoes an initial in- 
crease at the critical point after a microscopic time 

In our MC sweeps, the time evolution of M{t) is cal- 
culated through the definition 



M{t) 



1 

N 



I qMo-N 
\ 9-1 



(9) 



Here Mq = max(Mi, M2, • • • , M,) with Mi being the 
number of spins in the ?th state among q states. N = LP' 
is the number of spins on the suqare lattice. We use 
L up to 128. < • • • > denotes thermal averages over 
independent initial states and random number sequences, 
and [• • •] the disorder averages over quenched randomness 
distributions. The unit of time t is defined as a MC sweep 
over all spins on the lattice. 

The susceptibility plays an important role in the equi- 
librium. Its finite size behavior has been often used to 
determine the critical temperature and the critical ex- 
ponents ^/v and j9|. For the STD approach, the 
time-dependent susceptibility (the second moment of the 
magnetization) is also interesting and important. For the 
random-bond Potts model, it is defined as 



M(2)(i) 



1 

TV 



[(< M2(t) > - < M{t) >2)] . (10) 
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One studies its scaling behavior from mo = in short- 
time relaxation. Because the spatial correlation length 
in the beginning of the evolution is small compared 
with the lattice size i'' the second moment behaves as 
M(2)(t, L) ~ L-'^. Then the finite size scaling Eq.(|) in- 
duces a power-law behavior at the critical temperature, 



M^'^\t)^ty, y={d-f])/z. 



(11) 



In the above considerations the dynamic relaxation was 
assumed to start from disordered states with toq small 
enough. Another interesting and important process is 
the dynamic relaxation from a completely ordered state. 
The initial magnetization is taken exactly at its other 
fixed point mg = 1, where a critical scaling form 

M^^) (t, r, L) = 6" Wi^AfC^) {b-''t, b^/^r, b~^L) (12) 

is expected ||2^, |2^. This scaling form looks to be 
the same as the dynamic scaling form in the long-time 
regime, however, it is now assumed already valid in the 
macroscopic short-time regime. For the magnetization 
itself, b — i^/^ yields, for a sufficiently large lattice and 
at the critical point (t = 0), a power-law decay behavior 
of 



(13) 



The formula can be used to calculate the critical expo- 
nents ry = 2(3 /v or z. 

We further observe an evolution of the Binder cumu- 
lant U{t, L) = M^'^^{t, L)/M^{t, L)-l. A similar power- 
law behavior at criticality induced from the Eq.(|l^) 
shows that 



Uit,L) - t" 



d/z 



(14) 



on large enough lattices. Here, unlike the relaxation from 
the disordered state, the fluctuations caused by the ini- 
tial configurations are much smaller. In practical simula- 
tions, these measurements of the critical exponents and 
critical temperature are better in quality than those from 
the dynamical relaxation starting from disordered states. 



III. MC RESULTS 

In simulations, the heat-bath algorithm is used for MC 
updating, and up to 500 samples are chosen as disorder 
averages and 200-500 initial configurations and/or ran- 
dom number sequences are taken as thermal averages for 
each disorder realization (so total MC averages are over 
100,000). Statistical errors are simply estimated by per- 
forming three groups of averages with different random 
number sequences as well as independent initial states. 
We first focus our attention on evolutions of the magne- 
tization both from the initial states with small magneti- 
zation (rTiQ ^ 0) and complete ordered states (mo = 1). 

At the critical point (r = 0), achieved by generating 
the Boltzmann factors according to the distribution of 



FIG. 1: M{t) versus MC time plotted on a log-log scale and 
carried out on a 64^ lattice from the initial states mo ~ at 



and near the critical point, Kc and K, 



(±) 



(1.00 ± 0.01)A'c 




FIG. 2: U{t) versus MC time plotted on a log-log scale and 
carried out on a 64^ lattice from the initial state mo = 1 at 



and near the critical point, Kc and K, 



(±) 



(1.00±0.01)7^c 




Eq.(||), the curves M{t) should show power-law behavior 
as described by Eq. (^) or Eq. (|l^) dependent on the initial 
states. They, however, should deviate from the power-law 
scaling behavior for a small but nonzero t characterized 
by M(t,r) = mot^M(ti/''^T) or M{t,T) = t-'=^ F{t^/''''T) 
with respect to the mg ^ or mo = 1 initial states 
respectively. Now a question arises how to introduce a 
distribution which deviates from the critical point with 
distance r? We can first generate the Boltzmann factors, 
X = e~^, from the self-dual distribution Pxix) in Eq.(P, 



then modify them by the replacement x 



rl/(l+^) 



away 



from criticality |15| . We observe that the best power-law 
behavior appeared at the same critical point for both the 
random and ordered initial states. This fact indicates a 
second-order phase transition induced at as argued 
by Schiilke and Zheng |30) since Korder < Kc < Krandom 
if a first-order phase transition happened. We exhibit in 
Fig.l and Fig. 2 the double- logarithmic plots of magneti- 
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TABLE II: The values of indices ci, c„, 9 and y and estimated 
results z and r] estimated from the mo = 1 state for the g = 8 
on lattice sizes . 





mo = 1 


mo ~ 


results 


index 


Cl 




y 


z 


■n 


L 


= 32 


.0468(6) 


.603(7) 


.254(3) 


.597(6) 


3.32(5) 


.310(5) 


L 


= 64 


.0450(6) 


.597(7) 


.238(3) 


.526(6) 


3.35(5) 


.303(5) 


L 


= 96 


.0449(5) 


.593(7) 


.234(4) 


.508(5) 


3.37(5) 


.302(5) 


L 


= 128 


.0445(5) 


.586(7) 


.230(4) 


.497(5) 


3.41(6) 


.304(5) 



FIG. 3: Time evolution of magnetization for different lattice 
sizes from the ordered state at the critical point, plotted on a 
double-log scale. It shows the finite size effect on the M{t). 



zation versus the MC time on a i = 64 lattice to show 
this behavior obviously. 

Next, we carry out simulations at the exact critical 
point Kc to investigate the magnetization and Binder cu- 
mulant, starting from niQ — 1 ordered state. Such state 
is a ground state where all the spins orient in the same 
direction (e.g., cr^ = 1, z = 1,A^). In Figs. 3 and 4, M{t) 
and U{t) are presented and values of the scaling indices 
Cl = P/vz and c„ — d/z can be estimated from slopes 
of the curves in t =(100,1000). The results are listed in 
Table 2. Here we find that, unlike for relaxations from 
the disordered state, the fluctuations caused by ordered 
initial configurations are much smaller. As a result, these 
measurements of the critical exponents based on Eqs.(|l^) 
and (^ are better in quality than those from disordered 
states on Eqs.(|t) and (|ll|). So we take only 200 MC 
configurations for the thermal averages in the MC simu- 
lations. Now, from the results of ci and Cu, the critical 
exponents rj and z can be calculated based on Eqs.(p^) 
and (|lj), and their values are summarized in Table 2. 
To estimate the critical exponent 9, the random initial 
states are prepared with small magnetizations toq = 0.04 
- 0.01. The curves of M{t) for different mo on 128^ lat- 
tice with a double-logarithmic scale are displayed in Fig. 5 
which exhibits the finite tuq effect on M{t) and a very 
nice power-law behavior at toq — *■ 0. Thus the exponent 9 
can be estimated from slopes of the curves in t=(50,1000) 
by using Eq.(|^) after an extrapolation to mo — > 0. Its 
values are shown also in Table 2. 

As mentioned above, for a determination of the dy- 
namic exponent z and the magnetic exponent rj a dy- 
namic process starting from the ordered state (mo = 1) 
is more favorable since the thermal fluctuation is much 
smaller than for random states. Therefore we have 
mainly studied the short-time dynamic process from such 
an ordered state to calculate z and rj. On the other hand 
it is also desirable to measure these exponents from the 
mo = initial states p4[ . In our work we now take the 
evolution for the second moment of magnetization (t) 
from the mo = to verify the results of 77 and z. By 
input of values of rj and z already obtained from the in- 
dices Cl and Cu, we ckeck if they also satify the relation 
y = {d — rf)/ z in Eq. (pT|) . As expected we find they are 
consistent with less than 2% error for L > 96. 




FIG. 4: Time evolution of Binder cumulant for different lat- 
tice sizes from the ordered state at the critical point, plotted 
on a double-log scale. The tmic is extended to about one 
hundred. 




IV. CONCLUSIONS 

We have made a first attempt to apply the STD ap- 
proach to the 2D RBPF with a continuous distribution 
of quenched disorders. Our simulations verified that a 
second order phase transition has been induced. Then 
the universal dynamic scaling behavior in the short-time 
dynamics is used to estimate the exponents 9 and z, as 
well the magnetic critical exponent 77. It is found that 
they really violate the Ising-like universality class where 
9 = 0.191(1) and?? = 0.240(15) [|4|. Our resufl of the dy- 
namical exponent z = 3.41(6) is much bigger than that 
for systems without disorder, and it causes the micro- 
scopic time t„iic to be extended up to one hundred for 
the Binder cumulant from the ordered initial state. The 
value for the exponent 77 nearly completely overlaps with 
those listed in Table 1 , except that it is somewhat smaller 
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FIG. 5: Time evolution of magnetization for different initial 
magnetizations mo, plotted on a log- log scale. The curves of 
mo = O.OI and mo — 0.00 show nearly complete overlap. 



than that given in Ref . ||15| . More accurate calculation re- 
mains to be done to solve this controversy. 
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